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ABSTRACT 

X-ray observations unveiled various types of radio-silent Isolated Neutron Stars (INSs), phe- 
nomenologically very diverse, e.g. the ^Myr old X-ray Dim INS (XDINSs) and the ^kyr old 
magnetars. Although their phenomenology is much diverse, the similar periods (P=2-10 s) 
and magnetic fields (« 10^"* G) suggest that XDINSs are evolved magnetars, possibly born 
from similar populations of supermassive stars. One way to test this hypothesis is to identify 
their parental star clusters by extrapolating backward the neutron star velocity vector in the 
Galactic potential. By using the information on the age and space velocity of the XDINS RX 
J1856.5— 3754, we computed backwards its orbit in the Galactic potential and searched for its 
parental stellar cluster by means of a closest approach criterion. We found a very likely asso- 
ciation with the Upper Scorpius OB association, for a neutron star age of 0.42 ± 0.08 Myr, a 
radial velocity V^^ = 67 ± 13 km s^^, and a present-time parallactic distance djj?^ = 123l]^5 
pc. Our result confirms that the "true" neutron star age is much lower than the spin-down age 
ijsd =3.8 Myrs), and is in good agreement with the cooling age, as computed within stan- 
dard cooling scenarios. The mismatch between the spin-down and the dynamical/cooling age 
would require either an anomalously large breaking index (n ^ 20) or a decaying magnetic 
field with initial value Bq ss 10^"* G. Unfortunately, owing to the uncertainty on the age of 
the Upper Scorpius OB association and the masses of its members we cannot yet draw firm 
conclusions on the estimated mass of the RX J1856.5— 3754 progenitor 
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1 INTRODUCTION 

X-ray observations performed in the last two decades unveiled the 
existence of Isolated Neutron Stars (INSs) which are mostly radio- 
silent and not powered by the star rotation, such as the Myr-old X- 
ray Dim INS, or XDINSs (Turolla 2009), and the kyr-old magnetar 
candidates (Mereghetti 2008). Despite their different phenomenol- 
ogy, with XDINSs featuring stable and thermal soft X-ray emission 
{kT ~ 50-100 eV) and magnetars featuring transient and burst- 
ing high-energy activity and non-thermal spectral tails, both INS 
classes are thought to be linked by evolution. This is implied by 
their close locations in the period-period derivative (P-P) diagram 
and their similar values of the surface magnetic fields Bsurf, in- 
ferred either from the neutron star spin down or by the observa- 
tion of absorption features in the X-ray spectra, that can be both in 

* E-mail:rm2@mssl. ucl.ac.uk 



the ^ lO^-lO^" G range. While finding more similarities in their 
multi-wavelength emission would strengthen such a link, confirm- 
ing it is more challenging. 

If linked by evolution, both magnetars and XDINSs should 
follow a common evolutionary path. Interestingly, some magnetars 
seem to be associated with super-massive star clusters (Muno et 
al. 2005), thus with putative progenitors of > 4OM0, more mas- 
sive than those of "normal" neutron stars (8-25 Mq ; Heger et al. 
2003). However, identifying the parental clusters of the XDINSs is 
more complicated because of their larger ages with respect to the 
magnetars. Indeed, the results are affected by the uncertainty on 
the orbit extrapolation in the Galactic potential for time scales of 
a few Myrs. This depends on unknowns like the neutron star dis- 
tance and proper motion (hence its tangential velocity), measured 
from optical astrometry once the counterpart is known, and the neu- 
tron star age, inferred either from the spin-down (Tgd), once the pe- 
riod derivative is measured, or from the neutron star cooling (tcooi), 
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once both a model and a reference value of the surface temperature 
are assumed. A further unknown is the neutron star radial velocity, 
whose uncertainty dramatically increases the chances of spurious 
matches with clusters or stellar associations and hampers all iden- 
tification attempts (e.g., Mignani et al. 2010). However, in some 
cases the radial velocity can be inferred from the modelling of the 
bow-shock profile formed as the INS moves supersonically through 
the interstellar medium (ISM) and its fitted inclination angle with 
respect to the line of sight (LOS), like it had been done for Geminga 
(Pellizza et al. 2005).The association with the parental cluster also 
gives an estimate of the neutron star dynamical age (Tdyn)- 

The XDINS RX 11856.4-3754 (Walter et al. 1996) is the best 
target for this goal. It has the brightest optical counterpart amongst 
XDINSs {V ~ 25.7; Walter & Matthews 1997) and its proper 
motion and parallactic distance have been measured with high ac- 
curacy with the Hubble Space Telescope (HST), e.g. Walter et al. 
(2010), while the radial velocity has been inferred by modelling 
the bow-shock (van Kerkwijk & Kulkarni 2001a) detected in Ha 
by the Very Large Telescope (VLT). Moreover, the period deriva- 
tive of RX J1856.5 — 3754 has been measured (van Kerkwijk & 
Kaplan 2008) yielding the value of Tsd =3.8 Myrs. In this paper, 
we report on the search for the parental stellar association of RX 
J1856.5— 3754, based on the backward extrapolation of its orbit in 
the Galactic potential. The paper is divided as follows: the descrip- 
tion of the orbit simulation code and its application to to the neutron 
star and to candidate parental clusters/OB associations is given in 
Sectn. 2, the results are presented in Sectn. 3 and discussed in Sectn. 
4, respectively. 



2 ORBIT SIMULATION 

For the Galactic orbit simulation we used the code of Vande Putte & 
Cropper (2009), already successfully applied in Rauch et al. (2007) 
and in Vande Putte et al. (2009; 2010) and we refer to these publi- 
cations for further details. As discussed in Sectn. 1, in order to ex- 
trapolate the orbit of RX J1856.5-3754 in the Galactic potential 
and localise its putative birth place, accurate measurements of its 
proper motion p^^, parallactic distance d^^, and inclination angle 
i of the space velocity vector along the LOS are required. Both the 
proper motion and parallactic distance of RX J1856.5— 3754 have 
been repeatedly measured through optical astrometry techniques 
with the HST. As seen from Table [T] all measurements of the RX 
J1856.5— 3754 proper motion agree within the quoted uncertain- 
ties. For the parallactic distance, we assumed the most recent value 
of Walter et al. (2010), which confirms the earlier measurement of 
Walter & Lattimer (2002) and is consistent with that of Kaplan et 
al. (2002), but it is more accurate^ The inclination angle i along 
the LOS has been measured by van Kerkwijk & Kulkarni (2001a) 
by fitting the intensity profile of the bow-shock detected in Ha by 
the VLT. This is i = 60° ± 15°, which means that the neutron 
star would not move far from the plane of the sky. The alternative 
ionisation nebula model considered by van Kerkwijk & Kulkami 
(2001a), which yielded inclination angles closer to the LOS, has 



1 The first measurement of the parallactic distance, d^^ = 6llg (Walter 
2001), was not confirmed by Walter & Lattimer (2002) and Kaplan et al. 
(2002) and is not listed in Table 1, while those quoted in van Kerkwijk & 
Kaplan (2007) and Kaplan et al. (2007) were presented without any sup- 
ported evidence and are listed for completeness only. 



Table 1. Compilation of the proper motion (^^°) and parallactic distance 
(dJJS) measurements forRX J1856.5-3754. 



Parameter 


value 


Reference 




332 ± 1 


Walter (2001) 


(mas yr^l) 


333 ± 1 


Kaplan et al. (2002) 




331.2 ± 2.0 


Walter etal. (2010) 


(pc) 


140 ± 40 
117 ± 12 

167+15 


Kaplan et al. (2002) 

Walter & Lattimer (2002) 

van Kerkwijk & Kaplan (2007)* 

Kaplan et al. (2007)* 

Walter etal. (2010) 



This value has not been used in the current analysis 



been ruled out by Kaplan et al. (2002) because the nebula's open- 
ing angle would be incompatible with any of the parallactic dis- 
tance measurements. 

Using our orbit simulation code, we then calculated several 
test sets of neutron star tracks by looping on various values of the 
neutron star distance, the inclination angle (hence of the tangen- 
tial and radial velocity), and the neutron star age, computed around 
their reference values. In particular, we considered a grid of dis- 
tance values which are ±3cr around the best fit value of Walter et 
al. (2010), sampled with a spacing of 10 pc. From the correspond- 
ing sampled values of the transverse velocity, we then computed a 
grid of values for the radial velocity for different values of the incli- 
nation angle i (45°-75°) sampled with a 5° step. For the backward 
extrapolation time, we considered a grid of values for the neutron 
star age sampled at intervals of 10 kyrs. Since the uncertainty on 
the RX JI856.5— 3754 proper motion is < 0.6%, we neglected its 
influence on the orbit extrapolation and the neutron star birth place 
localisation. In the orbit computation, we did not account for pos- 
sible changes in the neutron star spatial velocity caused by close 
encounters with other stars or known star clusters. 

As a reference for the search of the parental stellar association, 
we firstly considered a sample of open clusters (OCs) selected from 
the latest version (3.1, released on 2010 November 24) of the data 
compilation of Dias et al. (2002), also referred to as the DAML cat- 
alogue. This contains entries for 2140 OCs, with sky coordinates, 
proper motion, radial velocity, and associated errors, together with 
information on the cluster metallicity, size, colour excess, Trumpler 
type, and age. Distances in the DAML catalogue are reported with- 
out associated errors, although a fiducial 10% uncertainty is proba- 
bly adequate for most cases (see, e.g. Vande Putte et al. 2010). We 
also used as a reference the catalogues of nearby OB associations 
of de Zeeuw et al. (1999) and Mel'nik & Dambis (2009), based on 
Hipparcos data. Both catalogues contain sky positions and proper 
motions in Galactic coordinates, radial velocities, and errors. In the 
de Zeeuw et al. (1999) catalogue, radial velocities are reported with 
no associated errors and we assumed a fiducial 20% uncertainty. 
Distances are derived from the trigonometric parallaxes of associa- 
tion star members (if available) or from photometric parallaxes. For 
the OB associations in the Mel'nik & Dambis (2009) catalogue 
we assumed their quoted 6% uncertainty on the distance, while 
for both proper motion and radial velocity (for which no error is 
given) we assumed a conservative 50% uncertainty. In all cases, we 
selected entries with non-null values of distance, proper motion, 
radial velocity (and associated errors). Although objects younger 
than ^ 100 Myr are obviously the most interesting candidates, ini- 
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tially we did not apply any selection based on the age of the OC 
or OB association, which can be uncertain by up to ~ 40%, as 
well as on other parameters, like the metalUcity or morphological 
type. Instead, we decided to use these parameters as validation ele- 
ments once a potential association was found. Our list include 439 
OCs from the DAML catalogue and 77 OB associations from the 
de Zeeuw et al. (1999) and Mel'nik & Dambis (2009) catalogues. 

We extrapolated back in time the orbits of the candidate 
parental OC and OB associations over the same age range as for RX 
J1856.5— 3754, using as a reference their nominal values of dis- 
tance, proper motion, and radial velocity. Then, we looked for the 
combination of neutron star parameters (age, distance, radial veloc- 
ity) which yielded the closest approach of the RX J1856. 5—3754 
orbit. We then regarded as likely associations those for which the 
approach was closer than a given threshold, defined as the overall 
tmcertainty on the computed separation. The association threshold 
accounts for the uncertainty on the cluster/OB association orbit ex- 
trapolation due to the random errors associated with their distance, 
proper motion, and radial velocity. As done in Vande Putte et al. 
(2010), we estimated this uncertainty through a Monte Carlo sim- 
ulation. For each object, we simulated 1000 different values of the 
distance, proper motion, and radial velocity, sampled within their 
formal errors, and computed the root mean square (rms) of the sep- 
aration between their backward-extrapolated positions at the refer- 
ence and the centre of the Galactocentric reference frame. In most 
cases, we found that the uncertainty on the orbit extrapolation was 
below 100 pc. The association threshold also accounts for the spa- 
tial extent of the OC or OB association computed from its angular 
size and distance and assuming, as a first approximation, a spherical 
symmetry. Since most OCs and OB associations tend to have irreg- 
ular morphologies, this is the most conservative assumption we can 
make. For simplicity, we did not account for two opposite effects 
which could influence the actual OC or OB association angular size 
in the past: its expansion due to intrinsic member star proper mo- 
tions and radial velocities, and its evaporation due to star escape 
from the local gravitational potential, which would yield angular 
sizes smaller and larger that those measured at the present epoch, 
respectively. 



3 RESULTS 

For completeness, we initially explored an age range of 2.8-3.8 
M3T around the RX J1856.5— 3754 spin-down age, with the caveat 
that this is an intrinsically uncertain age indicator since it depends 
on both the initial spin period of the neutron star and the value 
of the braking index n, which has not been measured yet for RX 
J1856.5— 3754. For the assvuned range of parameters we could not 
find a likely cluster association for RX J 1856. 5 — 3754. Although, it 
is possible that its parental OC/OB association is not a known one, 
owing to the relatively small distance (« 1.5-2.5 kpc) travelled in 
2.8^.8 Myr it is unlikely that it has not been discovered yet. While 
it is also possible that it has been filtered out in the sample selection 
(see Sectn. 2), the most likely conclusion is that the explored age 
range is not representative of the RX J 1856.5— 3754 age. 

This conclusions is confirmed by our measurement of its 
cooling age TcooI, computed using as a reference the most recent 
measurements of its surface temperature obtained from X-ray and 
optical-UV observations. It has been previously suggested that the 
surface temperature of RX J1856.5— 3754 is non-uniform (Pons 
et al. 2002; Braje & Romani 2002; Truinper et al. 2004), as in 
other cooling INSs. The discovery of X-ray pulsations at a period 



of ~ 7 s (Tiengo & Mereghetti 2007) supports this picture. Indeed, 
in a recent analysis of archival XMM-Newton observations of RX 
J1856.5— 3754, covering a time span of almost 10 years, Sartore 
et al. (2012) found that a two blackbody (BB) model is statisti- 
cally favoured with respect to a single BB in order to describe the 
0.15-1.2 keV spectrum. The resulting BB temperatures are = 
62.4 eV and Tc = 38.9 eV for the hot and cold components, while 
the corresponding BB radii are Rh = 4.7 (d/120 pc) km and Rc = 
11.8 (d/120 pc) km, respectively. When extrapolated to optical-UV 
wavelengths, the combined emission of the two BBs is consistent 
with the optical-UV fluxes obtained from //5r photometry (Kaplan 
et al. 2011), and further supports the two-BB picture. This gives a 
luminosity (at infinity) of 31.4 ^ logL (crgs^^) 32.5, where 
the uncertainty is computed accounting both for the unknown view- 
ing geometry and the uncertainty on the distance (Walter et al. 
2010). This range of Ivuninosity values is fully consistent with that 
of the single BB, 31.2 < log L (erg s"^) < 31.9, and those re- 
ported in the literature (see, e.g. Burwitz et al. 2003; van Kerkwijk 
& Kaplan 2007; Walter & Lattimer 2010) and lead to the same con- 
clusions on the source age. Assuming a minimal cooling scenario 
(Page et al. 2006, 2009), such luminosities imply cooling ages of « 
0.1-1 M5T, depending on the star mass and chemical composition. 
This is incompatible with the estimated spin-down age 4 Myr) 
of RX 11856.5—3754, unless the neutron star is closer than 90 pc, 
which is only marginally consistent with the 3a uncertainty on the 
parallactic distance (Walter et al. 2010). To svunmarise, the esti- 
mates of the cooling age Tcooi suggest neutron star ages as low as 
~ 0. 1 Myr Thus, wc reran our orbit simulations exploring the dy- 
namical age range 0.1-2.8 Myr. 

We found that the closest OC associations are with NGC 6475 
and ASCC 99 (see Table 2). For the former, the closest separation 
Ar — 160-170 pc for a neutron star age Tdy„ = 0.1-0.3 Myr, 
a present-time distance d^^ = 140-150 pc, and a radial velocity 
V^^ = 60-180 km s~^, while for the latter the closest separation 
Ar = 150-160 for Tdyn = 0.1-0.2 Myr, d^^ = 140-150 pc, and 

= 60-210 km s^^ In both cases, the values of the closest sep- 
arations are well above their corresponding association threshold 
(~ 15 pc). Moreover, the estimated OC ages of ~ 0.3 and ~ 0.5 
Gyr (Dias et al. 2002) would be much larger than the inferred neu- 
tron star dynamical age Tdyn of 0. 1-0.3 Myr and that of its putative 
massive progenitor (30-60 Myr). Thus, we ruled out the possible 
associations with NGC 6475 and ASCC 99. Similarly, we found 
possible associations with the Upper Scorpius, Upper Cen Lupus, 
and Lower Cen Crux OB associations. For the Upper Cen Lupus, 
the closest approach separation Ar ~ 55 pc, for Tdyn = 0.3-0.5 
Myr, d^^ = 110-150 pc, and V^^ = 50-180 km s"\ while for 
the Lower Cen Crux we derived Ar ~ 95 pc for Tdyn = 0.3-0.4 
Myr, d^^ = 90-110 pc, and V^^ = 90-190 km 8"^ However, 
only for the Upper Scorpius the closest approach distance is be- 
low the corresponding association threshold (~ 25 pc). We note 
that this value is dominated by the size of the OB association (de 
Zeeuw et al. 1999) rather than by the uncertainty on its orbit ex- 
trapolation, which is quite low thanks to the small uncertainties on 
the proper motion, distance, and radial velocity derived by Hippar- 
cos (see Table 2). In particular, for the Upper Scorpius we obtained 
Ar = 5-25 pc for Tdyn = 0.3-0.5 Myr, d^^ = 120-150 pc, and 
V^^ = 50-150 km s"'^. The inferred range of values for the neu- 
tron star present-time distance is consistent with that derived from 
the parallax measurement of Walter et al. (2010), d^^ = 1231^5 
pc. Conversely, fixing the neutron star present-time distance to this 
value would yield Tdyn = 0.42 ± 0.08 Myr, V^^ = 67 ± 13 km 
s"'^, and a closest approach separation Ar = 17 ± 8 pc, after ac- 
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Table 2. Name, coordinates, right ascension and declination proper motions (yUa; Ms), distances d, and radial velocities Vr for the OC and OB associations 
in the Dias et al. (2002) and de Zeeuw et al. (1999) catalogues which are potential birth places for RX J1856.5— 3754 (See Secto. 3). Values of the closest 
approach Ar and the corresponding neutron star dynamical age (Tdj,„), distance (d^®), and radial velocity (V^^) are given in the last foiff colvmms. 



Name 


RA 


Dec 


IJ-a 




d 


Vr 


Ar 










{hms) 


(C , „) 


(mas yr~i) 


(mas yr~i) 


(pc) 


(kms-i) 


(pc) 


(Myr) 


(pc) 


(kms-i) 


NGC6475 


17 53 51 


-34 47 36 


-1-1.67 ±0.20 


-3.60 ± 0.20 


301 


-15.53 ± 1.04 


160-170 


0.1-0.3 


140-150 


60-180 


ASCC99 


18 49 50 


-18 43 48 


+6.90 ± 0.64 


-2.50 ± 0.49 


280 


-31.29 ±0.40 


150-160 


0.1-0.2 


140-150 


60-210 


Up-Sco 


16 12 03 


-23 25 09 


+11.04±0.01 


-23.32±0.14 


145±2 


-4.6 ± 0.92 


5-25 


0.3-0.5 


120-150 


50-150 


Up-CenLupus 


15 08 12 


-43 45 06 


-i-21.30i0.35 


-23.13±0.14 


140±2 


-H4.9 ± 0.98 


55 


0.3-0.5 


110-150 


50-180 


Lo-CenCrux 


12 18 52 


-57 05 29 


-i-33.50±0.11 


-8.90±0.09 


118±2 


-H12.0 ± 2.4 


95 


0.3-0.4 


90-110 


90-190 



counting for the uncertainty on the parallactic distance. Thus, we 
deem the Upper Scorpius OB association as the putative parental 
cluster for RX J1856.5-3754. 



4 DISCUSSION AND CONCLUSIONS 

An association with the Upper Scorpius was originally proposed, 
e.g. also by Kaplan et al. (2002) and Walter & Lattimer (2002), 
assuming their own measurements of the RX J1856.5— 3754 par- 
allactic distance and leaving its radial velocity unconstrained. Their 
inferred dynamical ages for RX J1856.5— 3754 were Tdyn ~ 0-4 
Myr and ~ 0.5 Myr, respectively, close to the values inferred in the 
current work. Our conclusion is also in line with the more recent re- 
sults of Tetzlaff et al. (2010; 2011). In particular, for the same value 
of the neutron star distance d^^ = 123+^5 pc (Walter et al. 2010) 
as assumed by Tetzlaff et al. (201 1), the dynamical age that we de- 
rive for the neutron star (Tdyn = 0.42 ± 0.08 Myr) is fully consis- 
tent with their estimate (0.46 ± 0.05 Myr). The neutron star radial 
velocity that we infer (67 ± 13 km s^^), though, is larger than that 
derived by Tetzlaff et al. (2011), 6^20 km s^^, with any effect due 
the difference in the assvuned proper motion value being negligi- 
ble. The difference in the inferred radial velocity (still below 3cr) is 
probably ascribed to the different approach in selecting parameter 
values for the orbit simulation. Tetzlaff et al. (201 1) sampled a uni- 
form radial velocity distribution in the range -250-250 km s~^, to 
consider also directions much closer to the LOS, as predicted by the 
ionisation nebula model. Instead, by assuming the more likely bow- 
shock model (see footnote 2), we only sampled values in the range 
50-200 km s^^, according to the constraints on the bow-shock in- 
clination angle along the LOS. The difference in the inferred radial 
velocity can, thus, explain a factor of 2 difference in the closest ap- 
proach separation (for the same value of the neutron star distance). 
A similar result as ours has been obtained by Bignone et al. (in 
prep.) using a different orbit simulation code and a Monte Carlo 
approach to estimate the probability of a parental cluster associa- 
tion. They also find that Upper Scorpius is the only candidate with 
a non-negligible probability of being the RX J1856.5— 3754 par- 
ent system and obtain a value of Tdyn ~ 0.3 Myr, compatible with 
ours. To summarise, the agreement between these three results rules 
out that the proposed association between RX J1856.5— 3754 and 
the Upper Scorpius is affected by systematic effects and makes it 
statistically robust. 

Identifying the Upper Scorpius OB association as the birth- 
place of RX J1856.5— 3754 yields a robust value of the neutron 
star dynamical age Tdyn = 0.42 ± 0.08 Myr, which is almost 



an order of magnitude lower than the spin-down age Tsd = 3.8 
Myr. The dynamical age we derived is well compatible with the 
age, Tcooi, required to reproduce the observed source luminosity in 
standard cooling models (the minimal cooling scenario. Page et al. 
2006, 2009). "Fast" (or enhanced) cooling (e.g., Yakovlev, Leven- 
fish & Shibanov 1999 for a review) is not required to explain the 
observed properties of RX J1856.5— 3754. The mismatch between 
the spin-down and the dynamical/cooling age implies that the char- 
acteristic braking index of the source is n ~ 2Tsd/Tdyn ~ 20. 
Large and positive braking indices are routinely measured in radio- 
pulsars (e.g. Johnston & Galloway 1999). However, assuming that 
such a large n represents the average (constant in time) value in 
the standard expression for spin-down z> = —ku", where u is the 
star spin frequency and k a constant, would rule out any known 
form of braking torques. A possibility is to consider standard spin- 
down by magneto-dipole losses (n = 3 in the previous expression) 
but in the presence of a decaying magnetic field, which results in 
k OC decreasing in time. Self-consistent models for the coupled 
magnetic and thermal evolution of a neutron star were recently pre- 
sented by Pons, Miralles & Geppert (2009), to which we refer for 
details. Starting from a series of models with different mass M and 
initial surface dipolar field Bo, we looked for values of M and Bo 
which reproduce the present value of the luminosity L (two repre- 
sentative values 5 x lO^'^ and 10'^^ erg s^^ were chosen) at an age 
equal to the dynamical age, 0.42 Myr. This is achieved by models 
with Bo « lO^"* G and M in the range 1.1-1.7 Mq. Further im- 
posing that the period matches the observed one, P ~ 7 s (Tiengo 
& Mereghetti 2007) restricts the mass to be > 1.3 Mq. This further 
supports our conclusion that the true age of RX J1856.5— 3754 is 
sensibly shorter than the spin-down age and well compatible with 
our inferred dynamical age. The lower limit on the star mass when 
combined with the estimate of the star radius obtained by Sartore 
et al. (2012) from the two-blackbody fit, 12.5 < R (km) < 17.3, 
is rather non-constraining for the neutron star equation of state; we 
remark also that these estimates of M and R are model-dependent 
and should be taken with caution. A further possibility to reconcile 
the dynamical and spin-down ages is that RX J1856.5— 3754 was 
born with an uncommonly long period. Using the standard spin- 
down formula (P^ — Pg oc B^t), it easy to verify that to spin down 
the star at its present 7 s period in a time t = Tdyn ~ 0.42 Myr an 
initial period Pq ~ 6.5 s is required for B ~ 1.5 x 10^^ G. 

If, as it seems likely, RX J1856.5— 3754 originated ~ 
0.4 Myr ago in the Upper Scorpius OB association, we can con- 
strain the mass of the progenitor. The stellar population in Upper 
Scorpius has been extensively studied in the past (e.g. Priebisch et 
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al. 2002 and references therein). The estimated age of the associa- 
tion is ~ 5 Myr, according to mass estimates of Antares (an evolved 
Ml lb supergiant, ~ 20-25 Af©) and of the runaway OB star ( Oph 
(« 20Mo). Hence, the putative progenitor of RX J1856.5-3754, 
which went SN ~ 0.4 Myr ago, must have had a mass of 20- 
60 Mq, larger than usually expected for neutron star progenitors 
(8M0 < M < 25Mq - Heger et al. 2003) and close to those in- 
voked for magnetar progenitors (Muno et al. 2005). This, together 
with the high initial magnetic field required to reconcile the cool- 
ing and dynamical age in the framework of a decaying magnetic 
field, would then suggest that RX J1856. 5^3754 might have been 
an active magnetar in its youth. However, a recent re-assessment of 
the age of the Upper Scoipius association based on optical spec- 
troscopy (Pecaut et al. 2012) suggests a more likely age of ~ 11 
Myr. Hence, the masses of all members of the association would 
be lower than previously estimated (e.g. the mass of Antares would 
be f» I7M0) and the mass of the RX J1856.5-3754 progenitor 
would be ~ 18-20 Mq, i.e. still in the range expected for "nor- 
mal" neutron star progenitors. If one assumes that magnetars are 
indeed bom from supermassive progenitors, this might suggest that 
RX J1856.5— 3754 is not an evolved magnetar. However, any con- 
clusion has to be taken with due care. First of all, the estimates on 
the Upper Scoipius age and on the mass of its members are still 
debated, which makes it difficult to precisely constrain the mass of 
the RX J1856.5— 3754 progenitor. Second of all, any conclusion 
on the XDINS progenitors' masses must be verified against other 
XDINS/cluster associations. The only other XDINS for which both 
the optical parallax and proper motions have been measured (Ka- 
plan et al. 2007), and for which the Galactic orbit extrapolation can 
be reasonably well computed, is RX J0720.4— 3125. The source 
has been associated with the Trumpler 10 OB association (Kaplan 
et al. 2007; Tetzlaff et al. 201 1) for a dynamical age of 0.5-0.8 Myr. 
This would imply a progenitor no more massive than f» 2OM0, 
again close to that expected for "normal" neutron star progeni- 
tors. However, owing to the lack of information on the neutron 
star radial velocity from a yet undetected bow-shock, the associa- 
tion with Trumpler 10 is still tentative. The other way round, if one 
assumes that XDINSs are indeed evolved magnetars, the possible 
association of RX J1856.5— 3754 with a relatively low-mass pro- 
genitor might suggest that not all magnetars are bom from super- 
massive stars. Indeed, the proper motions of the magnetars XTE 
J1810-197 (Helfand et al. 2007), PSR J1550-5418 (Delleret al. 
2012), SGR 1806-20, and SGR 1900-^14 (Tendulkar et al. 2012), 
imply transverse velocities < 300 km s~^, which would not require 
hyper-energetic SN explosions from super-massive progenitors. 

To summarise, although it seems likely that the XDINS RX 
J1856.5— 3754 was bom ~ 0.4 Myr ago in the Upper Scorpius OB 
association, both the uncertainty on the age of the latter and the 
masses of its members make it difficult to draw firm conclusions 
on the actual mass of the neutron star progenitor, hence verify the 
possible evolutionary link between XDINSs and magnetars. Estab- 
lishing parental cluster associations for a larger sample of XIDNSs 
and magnetars, together with accurate studies of the cluster stellar 
population, is crucial to verify the connection between the proper- 
ties of their progenitors and the evolutionary tracks followed by the 
neutron star after the SN explosion. 
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